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ABSTRACT 

A  detailed  examination  of  the  literature  pertaining  to  subcooled 
critical  heat  flux  was  conducted.     Subcooling,  system  pressure,  mass 
flow  rate,  system  geometry  and  flow  stability  were  determined  to  be 
the  important  parameters  in  determining  critical  heat  flux.     The  effect 
of  each  of  these  is  discussed  in  detail. 

Various  schemes  for  correlating  experimental  critical  heat  flux 
were  subjected  to  critical  examination.     Form,  content  and  comparison 
with  experimental  data  were  the  bases  for  evaluation.     The  correlation 
of  Griffith  was  found  to  compare  most  closely  with  the  experimental 
data  over  a  broad  range  of  variables.     A  linearized  correlating  equation 
was  developed  for  the  MIT  low  pressure  (30  to  90  psia)  data,  including 
that  deterraned  as  part  of  this  study. 

Critical  heat  flux  was  determined  experimentally  for  pressures 
from  30  to  90  psia,  mass  flow  -ates  of  3.&9  x  10    to  lii.?6  x  10" 
lbm/ft^-hr  and  diameters  of  0.0li7  to  0.2h2  inches  over  the  subcooled 
region.     The  expe ~ime ntal  data  obtained  is  in  general  agreement  with 
that  of  other  investigations.     This  low  pressure  data,  when  combined 
with  the  results  of  previous  investigations,  provides  a  fairly  complete 
picture  of  subcooled  critical  heat  flux  for  tubes  up  to  0.25  inch 
diameter.     A  graphical  interpolation  scheme  is  presented  by  -which  it 
is  possible  to  obtain  predictions  of  critical  heat  flux  which  are 
generally  superior  to  that  of  the  mathematical  correlations. 

Thesis  Supervisor:    Arthur  E*  Bergles 
Assistant  Professor  of  Mechanical  Engineering 
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I.    IT. 

In  the  continuing  search  for  devices  of  higher  energy  flux  density, 
mere  and  more  emphasis  is  being  placed  upon  boili  a  means  of 

increasing  heat  transfer  rates  with  moderate  temperatur: 
Current  developments  in  high  energy  electromagnets,  high  performance 
electron  tubes  %nd  rocket  motors  require  very  high  heat  fluxes  in 
cooling  applications.     Although  heat  fluxes  well  in  excess  of  One 
Million  Btu  per  square  fcot  per  hour  can  be  obtained,  heat  fluxes 
cannot  be  increased  without  limit.     For  a  given  set  of  flow  conditions, 
there  is  a  limiting  Critical  Heat  Flux  (CHF),  at  which  the  transition 
from  nucleate  to  film  boiling  occurs.     At  this  point,  the  wall  teroerature 
increases  dramatically,  usually  resulting  in  catastrophic  thermo-mechanical 
failure  of  the  heat  transfer  surface.     One  also  finds  references  in  the 
literature  to  this  as  "burnout0,  and   "the  boiling  crisis". 

Often,  the  most  desirable  mode  of  boiling  for  heat  transfer 
applications  is  that  of  surface,  or  subcooled  boiling,  whe-e  the  average 
or  bulk  fluid  temperature  remains  below  that  of  saturation.     In  this 
case,  vapor  formed  at  the  surface  condenses  when  it  comes  in  contact 
with  the  cooler  main  stream.     The  Pressurized  Water  Reactor  uses  water 
at  high  pressure   (2000  psi)  as  coolant  and  moderator .     The  high  pressure 
allovs  sufficient  coolant  temperatures  for  further  heat  transfer  and 
power  generation,  while  maintaining  bulk  temperatures  below  saturation. 
The  iraintenance  of  subcooling  is  important  in  reducing  void  fraction 
and  retaining  continuity  of  the  moderator.     The  electromagnets  and 
electron  tubes  are  generally  cooled  by  low  pressure  systems  where  the 
expense  ind  complications  of  high  svstem  pressures  can  be  avoided. 
Utilization  of  subcooled  boiling  eliminates  the  complexities  of  two- 
phase  flow,  as  well  as  the  mechanical  complication  of  condensing 
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equipment;  however,  the  CHF  problem  is  still  p-esent  in  subcocled 
boiling , 

Over  the  p«  veral  years,  numerous  studies  have  been  concerned 

with  this  CHF  problem,     A  great  deal  of  experimental  data  is  available, 
under  various  conditions  of  flow.     Investigators  have  usually  selected 
certain  parftHetJtwa  and  Investigated  the  effect  of  v  these  within 

-tain  limits.     The  majority  oi  the  previous  experimentation  h 
been  conducted  *-ith  uniform  heat  flux  and  round  tubes.     This  investigation 
has  been  limited  to  this  region. 

The  first  phase  of  this  investigation  consisted  of  a  search  of  the 
literature  to  determine  the  flow  conditions  which  had  been  investigated. 
Techniques,  generally  in  the   form  of  correlating  equations,  have  been 
proposed  for  predicting  CHF,  by  almost  evervone  who  has  conducted 
inv  A  on  in  this  field.     One  phase  of  this  study  has  been  concerned 

with  a  critical  examination  of  the  techniques  and  their  form  and 
agreement  with  available  experimental  data. 

Another  phase  of  this  investigation  consisted  of  completion  oj 
some  previous  experimental  work,     Experiments  were  conducted  at 
selected  tube  diameters,  at  low  pressures  and  various  flow  rates. 
Data  thus  obtained,  when  combined  with  that  in  the  literature,  provides 
a  picture  of  CHF,  for  tubes  of  small  diameter  and  flow  rates  of  interest, 
over  th     complete  range  of  pressures,  from  30    >sia  tc  just  below  critical 
pressure. 
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This  was  a  two  part  stud*;;   s  of   the 

and  correlating  equations  and  experimental  determine 
Heat  Flux  (CHF)  under  certain  conditions. 

First,  citations  of  all  of  the  current  literature  concerned  with 
Subcode  d  Critical  Heat  Flux  were  obtained.     From  these,  the  ity 

of  the  literaoure  was  reviewed,  and  experimental  data  sed 

correlating  equations  were  recc  for  further  invc.  ion.     A 

large  part  of  the  work  in  this  field  has  been  accomplished  in  Russia, 
consequently  certain  information  is  available  only  in  the  Russian 
Journals,  untranslated.     It  was  usually  posrible  to  extract  data  and 
proposed  correlating  equations  without  complete  translation  of  the 
particular  paper, 

A  listing  of  proposed  correlating  equations  was  prepared.     These 
were  compared  with  actual  experimental  data  from  various  sources, 
including  the  experimental  phase  of  this  study,  in  an  attempt  to 
determine  the  validity  of  the  correlation  scheme.  . 

Experiments  were  conducted  to  determine  CHF  for  predetermined 
values  of  the  important  variables.     The  dependence  cf  C  on 

subcocling  was  determined  for  throe  fixed  values  of  rate, 

with  eras  cant  system  pressure,  tube  length  and  diameter*     This  data 
completes  the  low  pressure  end  oi  the  stue'v  of  Ornatski  a  higin 

(U,  !?  and  6),     Then,   using  one  oi  the  same  flow  rates  and  pressures, 
and  holding  exit  subcooling  constant,  dependence  of  CHF  upon  tube  diameter, 
and  upon  tube  length-diameter  ratio  was  determined,     Next,  the  variation 
of  CHF  with  changing  entrance  forms,  with  the  ether  variables  held 
constant,  was  investigated.     Finally,  tube  diameter  was  changed  to  obtain 
comparison  data,  and  to  determine  the  maximum  CHF  available  at  the 


, 


larger  diameter  with  the  experimental  apparatus. 

All  of  the  experimental  data  relating  tc  Subcooled  Critic; 
Flux  that  has  been  obtained  within  recenl  IT  was  collected, 

converted  to  similar  format  and  placed  on  machine  data  cards  in 
anticipation  of  future  studies. 

A  detailed  description  of  the  experimental  apparatus  and  prccec 
employed  is  contained  in  Appendix:  D  and  E. 
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III.    

FACT' __  K 

Although  there  has  been  a  ,  reat  deal  of  research  arid  irve  ion 

in  this  field,  particular?  in  the  1  cade,  there  exists  today  no 

r"1  larly  defined,  universal,  Mathematical  description  of  the  : 
which  results  ift  a  limiting  CKF.     Peel  Belling  is  itself  a  ccnplex 
problem,  the  ad  al  condition*  of  eubcooled,  forced  flow  ter 

to  cloud  the  picture  even  more*     Att  snpta  at  dimensional  analysis, 
and  treatment  of  the  problem  frcr  basic  equations  have,  in  the  r        , 
produced  little  information,  usually  yielding  long,  complicated 
functional  relationships  and  Involving  simplifying  assumptions 
necessarily  borne  out  by  experlr  nt. 

It  is  genersllv  accepted,  however,  that  subcocling, 
pressure,  mass  flew  rate,  flow  geometry  and  stability  are  the  m 

tors  in  determining  <"HF»     These  variables  evidence  various  decrees 
of  interdependence,  and  In  discus  be  effects  of  each,  this 

interdependence  must  be  borne  in     '      . 
Effect  of  Subccoling 

Subcoolinr  is  essentially  a  measure  of  the  aaottnt  of  pre-heating 
required  to  obtain  bi.  .     This  is  generally  measured  as   the 

difference  between  outlet  bulk  temperature,  or  more  precisely  enthalpy, 
and  that  at  saturated  conditions.     This  is  reasorable  :cr  uniform 
heating,  in  that  failure  normally  cr-urs  quite  near  the  downstream 
end  heated  section.     From  another  viewpoint,  cubccclin 

measure  of  the  ceding  effect  .     Thus  CH   increai    r 

with  subecoUng.     There  is  a  tendency  in  the 

as  linear  with  subccoling.     This  is  not  the  ea      ,   bat  it  is  in  al, 

not  a  grossly  inaccurate  approxima tion.     Tn  the  region  of  low  values  of 


subcooling  the  data  sh  -  imum  va! 

Clfi 

-  e  plaj  wrtant  role  :; 

the  effect  is  1c  arsatie  than  in  Fool  .     In  a  re<  .nt 

study,  it  is  stated  that  a  variation  of  e  fror 

tc  1000  psia  has  little  significant  effect  on  CBFi     this  in  tubes  r 
from  0*3  inch  to  0,29  Inch  in  diameter,  :  >  rn© 

out  by  the  other  exp© .  imental  data  (9).     A  pressure  increase  from 
600  psia  to  1000  psia  results  in  a  reduction  of  CHF  h  as  2'    , 

Further  increases  in  pressure  above  1000  psia  result  in  further 
redactions  in  Gl   .       t.  lower  pressures  there  is  disagreement  ever  the 
effect  of  pres-uro,     McAdams  (32)  and  Gunther  (33)  report  no  effect  up 
to  163  psia,  while  Mirshak,  et  al,   (31)  show  a  pressure  dependence  over 
the  entire  pressure  range  censidered,     Bergles  (ll),  reports  a  \% 
increase  in  CHF,  fcr  an  increase  in  pressure  frcn   20  to  85  psia* 
Chirkin  and  lukin  (2ii)  demonstrate  that  increasing  pressure,  in  the 
low  pressure  region,   increases  CHF,     Their  pressure  effect  is  Raich 
stronger  than  that  of  other  investigators,  however,   it  is  new  felt  that 
most  oi  their  data  suffered  fro®  instability  problems  which  could  have 
clouded  the  real  picture. 

The  experimental  data  dees  suggest*  however,    ihat  lor  fixed  geesset 
and  subcocling  under  conditions  cf  stable  flow,  there  is  scuie 
pressure  fcr  which  CKF  attains  a  jnaximura  value.     The  pressure  at  which 

r  tnBXimum  value  cf  CHF  occurs  is  net-  clearly  defined,  being  dependent 
upon  mass  flew  rate,  an  riablee  as  well.     It  is  safe 

to  say  thcu  h,  that  under  conditions  of  reasonable  lie  ,  ixiinum 

-c     rs  soinev;here  between  iiOO  and  600  psia.     Unfortunately  data  in  this 
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on  is  sea  i 

. 
tat*on  i  sube'ooling  is 

especially  a1 

t Doling  is 

3  suggests  that  at  pressure r  near 
high  values  of  subcooling  (Ahsat  greater  ':  *5JA°)i  ' 

somewhat  independen  subcooling.     Other  high  pressure  data  shows 

-;  sam  ■■■■!  trend  (1:,  5  .      "his  would,  c  Jra 

further  investigation  to  b«  throughly  substantiated. 
Effect  of  -ass  Flow  Pate 

ne rally,  in  the  gubcooled  region,  inc-r©  low  rate 

results  in  an  5  .     This  i  ::ly  visualized  since  the 

flow  sweeps  .ibble-    away  frc  neat  transfer  surface 

y,  repl—  atively  cold  water.     However*  the 

rate  o.-  reased  flow  rate  is  pressure  depc. 

At  high  press  the  effect  of  increased  b  i     as  less 

pronounced. 

—     e  ere  try 
The   r:<aicrHy  of  the  research  ?nd  st.  field  hae  be 

to  round  t  ular  flui<  .     f]  .  gross  e: 

B  hewev     ,  . 

In  al,  for    iven  flow  rate,  is  greatest 

pcular  tube  s  *  ?y  becc  re 

ple3E#     "Swirl  flow"  is  one  exception  tc  .     The  ratio  <  .ed 
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perimeter  to  total  yetted  perimeter  1  '■   i*1 

the  case  of  "swirl  flow"  (1?). 

Investigators  disagree  whether  th  where         has 

no  effect.  roroshchuck  and  Lantsman  in  a        udy  (9),  conducted 
tests  for  four  tube  diameters  over  mcs  of  the  high  pressure  ran 
(735  psia  to  25C&)  psia),  with  exit  conditions  from  hifb  s     ing 

(Ahsat*200  BTUA^n)  to  bulk  boiling  (xexit*0»2).  This  study  showed 
CHF  decreasing  consistently  with  increasing  channel  diameter*  Earlier 
investigations  show  decreasing  CHF  with  decreasing  diameter.  Both 
Chirkin  an  Iukin  (2li)  and  Kafenhaus  and  Bocharov  (26)  report  this. 
Again,  this  is  now  felt  to  have  been  in  a  large  part  dae  to  flow 
instabilities  in  their  systems.  DuPont  (8),  Jfirabak  (3)  and 
Zenkevich  (lit)  note  no  diameter  effect.  Bergles  (Ll)  notes  a  strong 
inverse  dependence  of  CHF  on  diameter  for  small  diameter  tubes,  but 
he  shews  the  Zenkevich  Prediction  (lU)  as  a  possible  asymptote  for 
diameters  greater  than  0.3  inches. 

V&th  flow  rates  high  enough  to  eliminate  the  effects  of  natural 
convection  within  the  passage,  channel  orientation  has  little  effect 
on  CRT. 

The  effect  of  heated  length  has  not  been  completely  determined,, 

In  general,  the  length,  or  length-diameter  ratio,  is  not  important  if 

the  length  is  long,  or  L/D  is  large.  The  popular  conception  is  that 

CHF  Is  determined  by  purely  local  phenomena,  thus  as  long  as  the  flow 

is  fully  developed  and  entrance  effects  play  no  part,  the  length  is 

not  significant.  But  the  question  of  how  long  is  long  enough  has  not 

been  satisfactorily  resolved,  partly  because  of  the  interaction  of  the 

variables.  In  general,  the  Russian  studies  consider  an  L/D  ratio  of 

10  to  be  large  enough  to  eliminate  the  length  effect.  This  conclusion 
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has  not  been  ocmpletly  verified.     "Berries   (11 V,   in  a  low  press  'Udy, 

shewed   .  vpendenre  existing  for  L/D  on  the  order  ■  .      "  croshr 

and  Fried  (18),  nt  data   for  tubi  .    >  inch  diameter  at  1^00  psia 

in  CHF   ^or  L/r  of  lh,  £C  ",  over  a  fairly 

braod  range  cf  exit  conditions.     nuPont  (3)  sugpests  cons  of 

ire  «grad  lent  down  "he  nested  tube.     The  terns 
is  of  course  dictated  by  the  First  Lav  of  Thermodynamics,  .his 

•lent  really  reflects  is  a  1  -diameter  consideration. 

In  general,  CHF  decreases  for  increasing  h/T)9  but  for  >ressures9 

the  dependence  of  CHF  upon  L/D  becoTnes  less  significant  for  values  of 
I/D  greater  than  30, 

ct  of  Flow  Stability 
This  refers  to  pulsations  in  fie-  -e,  resulting  from 

compressibility  effects  in  the  heated  sections.     If  there  is  an  expander 
or  accumulator  upsteam  of  the  heated  section,  there  is  a  distinct 
possibility  of  producing  a  flew  oscillation  which  will  periodically 
re'uee  the  flow  ne-?r  the  tube  <  ips  cause  a  premature 

failure  of  the  heat  transfer  surface.     Borgles  (11)  points  cut  that  there 
are  two  types  of  instability*     flow  oscillation  due  to  c empress ibility 
and  flow  excursion  due  tc  system  characteristics.     Flow  conditions  Under 
which  these  pulsations  do  effect  CHF  are  considered  unstable.     The  result 
of  unstable  flow  is  a  drastic  reduction  in  rHF5  with  unstable  flow' 

heat  fluxes  as  low  as  20$  of  the  stable  value.     Stoble  flow 
usually  can  be  insured  for  subcode4  belli  tailing  some  sort  of 

thrott"  tartee  (orifice,  valve,  etc.)  in  t'  ten  just  up-stream  of 

I  section.     Poroshehuek  ran   ase  -ure  drop  of 

li70  psi  in  their  high  pressure  studies  (9,  10).     Perries  (ll)  reported  a 
pressure  drop  of  approximate  ll  in  his  low  pressure 
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Further  discussion)       stability  pr- 
Itef,  7,  11     •  .       '  ility 

system  re  .ubes  from  cc 

inlet  and  exit  headers,  where;  a  suba 

volume  may  exist.  In       t,         I  on   of  icing 

Jevie©|     -  »*resuliant  losa    osmstream  sysl         ,    -ore 
than  justified  by  the  increase 
Effect  of  Secondary  Factors 

Several  other  factors  have  been  investigated  to  determine  how  they 
effect 

rface  roughness ,  for  reasonably  smooth  (<       I  quality) 
onels  has  little  or  no   eiiect  (7,  13 )• 

Dissolved  gas  in  the  flow  i     ,  in  small  amounts,  has  a  n 
•■-■&  on  CRF  (17)* 

For  heat  flux  unevenly  distributed  areund  the  heated  perimeter, 
peak  heat  flux  was  found  to  increase.  For  heat  flux  unevenly  distributed 
alonp  the  length  oi        ,ed  channel,  qualitatively  similar  results 
have  been  obtaiood  (7),  Tests  have  been  conducted  under  renditions  where 
the  I     lux  -was  forced  into  a  cosine  distribution  along  the  tube  and 
CHF  was  shown  to  be  reduced  by  seme  30$  (11), 
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General  discussion 

TJoaierous  schemes  have  been  nror>osed  as  cor: 

uich  CFf  can  supposedly  be  predicted.     The  earl:5 

of  a  inctioi 

*  (qA)cr  -  A*  Be  Tf 

•where  B,   D  and  F  are  the  irnportant  variables,   or  rations  of  tl 

and  g  are  experiment ally  de  orients,       sually 

correlations  of  this  tyr>e  tei  tier  between 

variables  j  her  cases  the  exponent  tactions  of 

the  ccntr-  asters.     Another  brei  cry  includes  the 

"sjp-r  on"  correlations |  ing  ?c  lata  -with 

forces?  convection  heat  transfer  infc  .  are 

often  lonr  and  difficult  tc  use, 

Several  attempts  have  been  Trade  at  s  tt,  Venice vich 

(lii)  and  Fyabov  and  Bersina   (2  1  analysis  and  Griffith 

(13,  17)  with  a  theoretical  analysis.     Other  equr  -.weloped 

by  *bnite  frrce"  c  sis,  .        iiPont 

f avers  t   form: 

(q/A)cr  *  C  (I  ♦  OlJL)  (1  ♦  r2.3)   (l  =  fs. 

yhere  A,  B  and  D  are  the  influencing  les,   (pvr  ling, 

diar^eter,  velocity,)  and  C,  Cp   Gg,  etc.   are  constants. 

In  p   r-ral,  the  corrc  es  of  the 

ables,  departing  radically  fren  tl  "  in  the  excluded 

of 

an  experifflental  study,   using  cert  .       '  is 

•;ue  is  usually  aceompani  so 

manv  percent  1  results  agrees  within  plus  uch. 

1* 
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In  cases  such  as  this,  the  potential  user  is  obliged  to  examine  the 
filteri  redure  as  closely  as  be  dees  the  correlation  itself. 

As  a   result  of  this  state  of  affairs,  one  ph; 
was  concerned  with  examination  of  correlating  equations.     The  listing 
of  correlations  is  not  all  inclusive.     Several  correlations  were  examined 
and  discarded,  generally  due  to  their  naikotf  similarity  tc  otfaei    . 

Of  the  correlation  schemes  examined,  the  following  are  those 
which  agreed  reasonably  well  with  the  experimental  data,  or  illustrate 
some  point. 

Details  of  Particular  Equations 
Bellt   (reft  21) 

(q/A)cr  -  0.1,27  x  yfi  ,  Ja.59  .  o^1'81  *  (i^)  M»  x^.m 

Range  of  values:     (qA)CP  •  Btu/hr-ft2;  p  *  2000  psiaj  G  •  .2  x  10°  - 

5.0  x  106  lbm/hr-ft2|   Ahsat  -  0  -  130  Btu/lbm 

This  correlation  is  compared  with  experimental  data  (figure  9) 
for  conditions  very  nearly  identical,  from  three  different  sources. 
In  obtaining  this  data  from  the  various  references,  a  certain  variation 
was  unavoidable}  two  people  do  not  conduct  tests  at  precisely  the  sam« 
conditions.     However,  the  variation  was  kept  to  less  than  the  expected 
experimental  error.     The  data  from  reference  15  was  obtained  in  1951; , 
prior  to  the  general  acceptance  of  the  need  for  stability,  and  may  be 
of  questionable  value. 

It  can  be  seen  that  the  correlation  follows  the  general  tiend  of  the 
data,  with  variations  of  15$  to  30%. 
Jena  6  Lottes:  (ref:  1£) 

This  correlation  was  based  on  data  taken  at  UCLA  and  Purdue  over  a 

range  of  pressures  from  500  -  3000  psia.     The  values  of  K  and  n  were 
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determined  on  the  basis  of  the  data  such  to  make  all  deviations  positive. 
(q/A)cr  -  K  (0/l06)n  '(at^   *22 

UCLA  Data  correlated  t  -OS*,  +23$,  D  ■  .26  in. 

p     K     n 

$00   .817   .16 

1000    .626   .275 

2000    .Ui5   .50 

Purdue  Data  correlated  to  -0j6f  ♦60*.,  T1  ■  .3lt3  in. 

P  K  n 

3000        .915      .275 

2000       .5U5     .50 

3000         .30         .725 
The  correlation  was  plotted  for  data  (figure  9)  of  approximately 
D  •  .160  in.,  G  •  2.1  x  10^  lbm/hr-ft2,  at  p  -  2000  psia  usinr  the 
constants  and  exponents  of  the  Purdue  data. 
frrabov  and  3orzina t     (reft  36) 

Range  of  valuees     D  •  2  «■  •  .079  in.,  AT-a*  •  18  -  590  deg  1 , 

I/D  ■  20,  p  •  10  -  200  atmospheres  »  11(7  -  291:0  psia,  0  -  3.69  x  10    - 

22.11;  x  106  ltw/niH-ft2 

2)  (q/A)cr  -  (q/A)ph  *  2S7  V6AW6?  K  W^O*  l  -f>v)*515     I*1 


,(>r  t).jj  v^  ^ — 


Range  of  values:     D^lran  (.158  in.),       Tsa^  •  IS  -  k30  deg  F, 
I/D  -  10  -  220,  V  -  3  -  U8  ft/sec,  p  -  70  -  210  atm  -  1030  -  3090  psia 
oth  of  the  above  correlating  equations  were  developed  from 
dimensionless  equations.     In  each  case  (q/A)  .   is  the  same  and  is  obtained 
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from  a  plot  of  pool  boiling  CHF  derived  as  a  function  of  pressure 
correlating  equations  were  derived  from  a  dimensionless  equation; 
1)  K^  -  3710  IV6  KAT-67  Pr"#i*       >  -  2am  »  .079  in.) 
or  2)  KAq  •  2370  Y^*6  KaT»6?  ?r~^     (D  *  Iwm  *  .l£l   in.) 

where  KAq  *  &qA)cr  kfgpv  J  go 

■■>    '■■■•  ■  ■--  ' ' '  '  ■ 

ks  Tsat  ^1  -<°v)g 


Kw-V//9! 


^LEIzIjl. 


c-* 


K  T  ,  ATsat  hfg  ^»v  .rgp 


J 

g*-Tsat 


JEHfi. 


It  is  assumed  that  (q/A)cr  ■  (o/A)pb    *  (  Aq/A)er 

The  2  mm  (.079  in.)  correlation  was  derived  from  the  data  of  Omatski 
and  Kichigin  (U,  5,  &  6)  by  the  usual  technique  of  letting  the  data 
determine  the  coefficients  and  exponents.     Generally  the  data  falls  within 
-  30$  of  the  correlating  line.     The  correlation  was  also  confuted  for  the 
2  ma  data  taken  during  the  experimental  part  of  this  study.     In  this  case 
the  correlation  did  not  follow  the  basic  trends  of  the  data  and  was  in 
error  by  as  much  as  100$  at  some  subcoolings. 

The  correlation  for  D  *  lj  mm  (.158  in.)  *jas  computed  and  plotted  for 
G  -  2.1  x  10    lbm/hr-ft2,  p  ■  2000  psia  and  compared  with  the  data  for 
Poroshehuk  and  Lantsman's  (9)  k  mm  and  8  mm  tubes.     The  correlation  fit 
the  8  mm  (.316  in.)  data  very  closely  ,  but  was  K$%  -  30£  low  on  the 
k  mm  data.     It  was  not  plotted  on  figure  9  to  avoid  confusion  as  it 
fell  generally  on  top  of  the  correlation  of  Zenkevich  and  Subbotin. 
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Zenkevich:     (ref:  Hi) 

(q/A)cr  »  hfg  te<rs  gp\^    *  p.5  ♦  !8Mhsat1     x  10~$ 

Range  of  values:     p  »  Hi  -  165  psia,    Ahgat  *  20  -  2        Etu/lbm, 

G  -  ,22  x  106  -  9.96  x  106  lbm/hr-ft2,  D  «  .16  -  .50  in. 

This  correlation  is  linear  with  subcooling  as  comonly  noted  for 
Rusaiin  analyses.     Specific  instructions  as  to  the  temperature  to 
evaluate  viscosity  at  are  lacking.     Using  the  fluid  saturation  temp- 
erature yields  an  equation  linear  in  subcooling,  while  using  fluid  bulk 
temperature  produces  a  slight  curvature  in  the  resulting  plot.     This 
correlation  is  plotted  in  figure  8  for  both  conditions  and  it  can  be  seen 
that  the  equation  does  not  really  reflect  the  experinsental  results. 
Zenkevich  and  Subbotin:      (ref:  lit) 

(qA)cr  -  397  GiATsat«33  (Tg/Vfg)-1.8      ^^.^2 

Range  of  values:    ATsat  *  18  -  180  deg  F,  B  »  U  -  12  mm  *  .158  -  .U72  in. 

p  ■  mo  -  210  atmospheres  -  2060  -  3090  psia. 

This  correlation  was  plotted  for  G  ■  2.1  x  10^  lbm/hr-ft2  for  2000 
psia  and  2500  psia  and  compared  with  the  data  of  Doroshchuck  and  Lantsman 
(9)  for  tubes  of  h  ram  (.158  in.)  and  8  ram  (.316  in.)*     At  2000  psia  the 
correlation  does  an  almost  remarkable  job  of  predicting  CKF  for  D  *  8mm. 
At  2^00  psia  the  correlation  is  very  close  to  the  experimental  8  mm  data 
although  approximately  20$  low  at  low  subcoolings.     The  correlation 
does  not  fit  the  k  mm  (.158  in.)  data  as  well,  with  differences  of  1$%  - 
30$  at  both  pressures.     Figure  9  shows  the  correlation  at  2000  psia  for 
the  h  mm  diameter. 
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Gambill:     (ref s  30) 

.  r— » 

►923 


(q/Mcr  *  !   hfgfv/^go  a  (f  1  -I°v) 


# 


A 


♦  K«  (k/D)  Ren  Pr11  (t„  -  tb>Bo 
where  K  •  #12  -  .17*  K1,  m,  &  n  depend  on  the  non-boiling  turbulent 
flow  correlation  used. 

Range  of  values:     V  *  .05  -  17k  ft/sec,  p  •  2u2  -  3000  psia, 
ATsat  -  0  -  £06  deg  F,  a  *  1  -  57*000  g»s 

This  correlation  is  of  the  superposition  form  where  a  subcooled 
pool  boiling  burnout  term  is  added  to  a  non-boiling  forced  convection 
term |  typically  the  well  known  Kc  Adams  correlation  for  turbulent  flow. 

Gambill  leaves  some  latitude  for  the  evaluation  of  the  constant 
eoefiiecients  and  exponents.  Evaluation  of  the  wall  temperature  (ty) 
in  the  forced  convection  term  is  done  by  an  auxiliary  plot  of  >^Tsa^ 
vs#  Tsat/^crit  based  on  Bernath's  correlation  of  wall  temperature  at 
burnout  (16).  This  is  explained  in  detail  in  the  section  on  Bemath's 
correlation.  The  correlation  was  plotted  under  different  conditions 
with  the  following  results: 

1)  D  -  .2h2  in.,  G  *  3.69  x  106  lbm/ft2hr,  p  *  90  psia  (figure  8) 
The  correlation  fell  about  35$  lower  than  the  experimental  data. 

2)  D  •  .079  in.,  G  »  3.69  x  106  lbm/hr-ft2,  7.38  x  106  lbm/hr-ft2, 
&  lb.76  x  10    lbm/hr-ft2,  p  ■  90  psia  (figure  l*)s     The  correlation 
fell  as  much  as  100''  lower  than  the  experimental  data. 

3)  P  -  .180  in.,  p  «  2000  psia,  G  »  2.1  x  ID6  lbm/hr-ft2  (figure  9) 
The  correlation  is  conservative   'low),  but  does  not  show  the  trends 
of  the  data  completely,  nnd  is  in  error  by  as  much  as  2£  '• 
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Berrath:     (reft  16) 

'i      i    ■■  »■' 

(q/A)CP  -  hbo  ^*w  "  *bH>o 

where:    h^  »  5710  +  U8  ¥  D"*6      Btu/hr-.ft2«deg 

and:     t^  -    32  *  (9/5)  f  57  ln(p)  -  Skip)     -  Y  j     deg  F 

Range  of  values:     Round  tubes,  p  •  $00  •  3000  psia,  V  »  i*.     - 

5h  ft/sec,  where  D  is  measured  in  feet,  and  p  in  psia. 

This  correlation  is  based  on  a  limiting  heat  transfer  coefficient 
and  a  correlation  of  wall  temperatures  at  burnout.     The  correlation 
was  plotted  (figure  9)forG-2.1x  106  lbm/hr-ft2  at  2000  psia,  and 
D  ■»  .180  in.  for  comparison  with  experimental  data.     The  correlation  is 
slightly  conservative  predicting  this  data  within  approximately  $%9 
tist     (reft  17) 

(<l/A)cr  -  106  x  0.28(  h^lO3  )"2^  (  1  +  0/10*7  )2  9-.<*n2J/V    Btu/hr-ft2 

Range  of  values:     G  ■  0.2  x  106  -  8*0  x  106  Ibm/hr-ft2,  p  *  l8£0  - 

21^0  psia,    Ahgat  *  0  -  180  Btu/lbm,  l/D  •  21  -  365. 

This  correlation  is  compared  with  experimental  data  on  figure  9. 
It  can  be  seen  that  it  follows  the  general  trend  of  experimental  data, 
with  variations  of  15$  to  1*0$. 
Griffith:     (ref:  13,  17) 

(q/A)cr  «  f(  p/pc  )   (F)     g  (hg  -  1%)  fg  (Z3!-^  )(   *i\ 

(F)  =  1  ♦  (Jfcfc  x  ID"6  )  -  0.1k  Kj_  ♦  0.5  x  10"3  (Re^  Kx  )^ 

and,  f(  p/pc)  is  presented  graphically  in  the  reference. 
Range  of  values:     p  ■  lii.7  -  3000  psia,  7  ■  0  -  110  ft/etc, 
^Tsat  -  0  -  280  deg  F, 
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The  correlation  was  plotted  and  compared  with  data  at  two  res: 

90  psia  (figure  8)  and  2000  psia  (figure  9).     At  90  psia  the  corral 
is  outstanding  and  does  an  excellent  job  of  predicting  the  data. 
2000  psia  the  correlation  is  higher  than  the  data,  yet  only  by 
■'irshak,  et  al:     (ref:  13) 

(q/A)cr  =  h §0,000  (   1  ♦  .0131  p  )(  1  ♦  .00503  ATsat)(l  +  .036$  V) 
Range  of  values:    p  *  2£  -  8£  psia,  ¥  *  $  -  h$  ft/sec,  D  -  .21  - 
♦h6  in.,  whore  £Tsat  is  measured  in  dog  F,  and  (q/A)cr  in  Btj^/hr-ft2. 
This  correlation  was  compared  with  experimental  data  for  .2U2  in. 
diameter  tubes  at  90  psia  (figure  0)  and  .30$  in.  diameter  tubes  at  30 
psia  (26).     In  both  eases  the  correlation  curve  was  conservative  (lower 
than  actual  data). 
DuPont:     (ref:  __ 

(q/A)cp  »  595,000  (1  +  ,0liV)(l  *  .0055  ^T^tHl  ♦  .»       ^  3Cg) 

Range  of  values:     V  -  8,$  -  1*0  ft/sec,  AT^  *  kO  •  120  deg  F, 

D  ■  .3  -  .68  in.,  p  e  60  -  1000  psia,   ATscg  ■»  Subcooling  Gradient  * 

(Texit  *  *inlet)  I**  2  foot  length  of  tube. 

DuPont' s  correlation  is  similar  to  the  one  developed  earlier  by 
Kirshak,  et  al  (13)j  a  combination  of  linear  factors  of  velocity  and 
subcooling.  In  the  DuPont  tests  a  two  foot  length  of  tube  was  used. 
Subcooling  Gradient  (4  ^gcg)*  the  temperature  difference  through  the 
tube,  was  the  last  factor  used  to  correlate  experimental  data,  using 
the  form  of  a  previous  correlation  for  annular  geometry. 
Doroshchuck  and  Lantsirgtn:(ref :  9) 

(q/0cr  -  (q/A)0  *  \(p)  [(.3l6/0)*33  .  l]     Btu/hr-ft2 

where:     (q/A)0  is  determined  from  a  plot  of  (q/A)cr  vs.  A^SSt^  fer 

8  ram  (.316  in.)  tubes  for  the  pressure  in  question. 

Range  of  values:     D  ■  .117  in.  -  .316  in.,  p  *  50  -  170  aten  •  73'3  - 

2$00  psia.     A(p)  is  a  tabulated  function  of  pressure. 
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p  -  atm 

pada 

738 

A(p)  -  Btu/hr-fiT 

50 

2.55  x  106 

80 

111*0 

2.1*1  x  106 

100 

11*70 

2.31  x 

ifco  % 

2060 

2.13  x  106 

170 

2500 

1.99  x  10 

:3ased  en  the  8  ram  (.316  in.)  data  of  Doroshchuck  and  Lantsman,  the 
correlation  predicts  critical  heat  flux  for  the  other  diameters  with  a 
maximum  deviation  of  30$,  for  fixed  values  of  the  other  variables.     "1th 
other  experimental  data  the  variation  is  slightly  greater. 
KIT  Small  hamster: 

(q/A)cr  ■  1.736  x  106  (1*.1*8  *  .Ol^h^M^)^'  (^r^ft^p12 
where:     GQ  -  106  lbm/hr-ft2 

D0  ■  0.01  in. 
Ranee  of  variables:     D  -  .05  -  .305  in.,  p  *  30  -  90  pala,  G"2x  10° 
15  x  106  Ibn/hr-ft2,  I/D  •  5  -  60,    A^t  *  0  -  150  Btu/lbm 
This  equation  was  developed  by  the  authors  for  the  low  pressure 
data  fron  MIT  and  fits    he  majority  of  the  data  to  1  2$%» 
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Scope    of  Experimental    Study 

As  mentioned,  there  is  a  grer»t  deal  of  experimental  data  available 
in  the  literature.     Although  each  investigator  has  attempted  tc  be 
somewhat  systematic  in  his  own  method  of  attack,  there  are  very  few 
cases  where  data^from  cne  investigation  can  be  compared  with  that  from 
another.     Generally  there  are  substantial  differences  in  the  variables, 
and  lacking  a  gcod  correlating  equation,  cor-parison,  except  in  a  very 
general  ;:av,   is  irpossible.     The  experimental  phase  of  this  study  was 
conceived  with  this  situation  in  rrind.     Ornatski  and  Kichigin  (h$  5 
and  6)  have  conducted  tests  with  tubes  of  two  millimeter  diarseter,  and 
four  mass  flow  rates,   over  b  range  of  pressures  from  li*?  psia  tc  28UO 
psia.     Part  of  the  present  study  was  to  complete  the  lew  pressure  end  of 
the  2  mm  map.     Also,  using  the  same  size  tube  and  flow  rate,  other 
-parameters  were  varied.     The  results  of  this  investigation,  when  combined 
with  the  data  from  Ornatski  and  Kichigin,  present  a  fairly  coirplete 
picture  of  CHF   for  tubes  of  two  millimeter  diameter. 

Experiments  were  conducted  to  yield  the  following  sets  of  data: 

a.  CHF  at  various  subcooling  for  three  mass  flow 
rates  (G  *  3.69  x  106,  7.37  x  106,  lh.76  x  106 
lbm/hr-ft     -  identical  tc  those  used  by  Ornatski 
and  Kichigin)  with  D  *  2™  (.079  in.)  for  exit 
pressures  of  30  psia  and  90  psia. 

b.  CEF  at  various  tube  diameters,  with  diameters 
varying  from  0.01*7  in.   to  0.2L;2  in.,  with  mass 
flow  rate  held  constant  at  3.69  x  10^  lbm/hr-ft2, 
and  at  constant  exit  subcooling  of  35  Btu/lbm. 
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i.     With  L/D  held  constant  at  a  value 
of  hO, 
ii.     With  L  lie  Id  constant  at  a  value 
of  3.15  inches 

c.  CIH    for  various  L/D,  for  the  same  mass  flow 
rate  and  exit  subcocling  as  in  part     b.   above 
with  D  •  0.079  in.  l/V  from  3.5  to  60. 

d.  CHF  at  various  subcoolings  for  the  same 
mass  flow  rate,  with  L/D  of  Uo  and 

D  *  0.2li2  in.     This  dat3  was  obtained  in 
order  to  present  a  comparison  with  the 
corresponding  D  *  0.079  data,  and  as 
substantiation  of  the  CHF  verses  diarreter 
data. 

e.  (TIF,  for  the  same  flow  rate,  exit  subcooling 
and  L/D,   D  »  0.09U  in.  with  changing  entrance 
conditions.     In  addition  to  the  calming 
length  form  used  in  the  majority  of  the 
study,  four  entrance  forms  were  investigated. 
In  all  cases,   the  entrance  was  placed 
immediately  up-str3am  of  tie  heated  length. 
Entrance  forms  investigated  were: 

i.     Sharp  Edge  Entrance. 
ii.     i4$«<?.0gr:  e  Chamfer. 
iii.     Rounded  Entrance,  radius  of 

entrance  equal  to  2  diameters. 

iv.     Reentrant,  length  of  reentrant 

section  equal  to  10  diameters. 
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f,     T'axirum  Heat  Flux.     An 

of  the  system  capabilities,  using 
0.212  in.      .    .  ,  with  l/V  of 

7,?£  an  ^  Liable  flow 

cit  subcooling,  tc  dete~ 
tn*fe  naximum  obt  obis 

diamet  l . 
Results 

A  total  of  117  CHF  data  points  "were  obtained  in  the  experimental 
phase  of  this  study.     Of  these,  twelve  were  obtained  using  the  Beryllium 
Copper  tubing,     A  detailed  description  c:  the  problems  encounter 
with  this  tub  inn  is  contain  endix  D. 

Figures  3  and  h  show  the  CHF  data  obtained  with  D  »  0.079  in,, 
,i-  the   three  mass  flow   rates.     These  figures  show  the  same  general 
#ndeacs  of  OH    on  mass  flow  and  subcooling  ss  the  majority  of  the 
other  experimental  data.     This  lex?  pressure  CHF  data  show  the 

-inimum  CHF   in  the  subecoled   region  more  clearly  than  does  the  data  at 
higher  pressures.     It  can  be  seen  that  imum  CHF  is  encount  at 

approximately  30  3tu/lbm  subccoling,  and  that  the  subecciing  for 
•  inimum  CHF  increases  for  increasing  mass  flow  rate,     CHF  beei  Me« 
nearly  indapendant  of  subcooling  at  the  lowest  flow  rate,  at  30  psia. 

The  dependence  of  CHF  upon  tube  diameter  is  shown  in  figure  5, 
Investigations  were  conducted  with  L/D  held  constant  and  vith  L  held 
constant.     The  crossing  of  the  two  curves  merely  shows  the  diameter  at 
which  the  selected  values  of  L  and  L/T)  coincide.     The  data  obtained 
with  L  held  constant  shows  a  greater  CHF  fcr  the  larger  diameters, 
reaching  a  maximum  increment  of  approximately  1$%  at  the  0.21+2  in, 
diameter  tube.     The  value  of  CHF  can  be  seen  to  become  less  dependent 
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upon  diameter  for  the  larger  tubes,     Limits  of  the  apparatus  precluded 
investigation  of  diameters  larger  than  0,2lj2  in.  at  the  flow  rate 
selecte  .his  study. 

The  comparison  between  CHF  data  at  D  ■  0,079  in,  an i  that  at 
D  *  0.2ii2  in,  can  be  seen  in  figure  6#     This  data  provides  an  evaluation 
of  the  validity  «f  the  curve  of  CU'  verses  diameter  (f  if  ore  b)  which 
was  determined  lor  a  particular  value  of  exit  subcooling.     It  can  be 
seen  that  proportional  difference  i:  xjtween  the  two  dia 

is  nearly  constant  over  the  range  cf  exit  subcoolinga  considered, 

Inv  ition  of  the  effect  cf  changing  length-diameter    i 

upon  CHF   (figure  7)  shows  that  Ci27   increases  for  decreasing  L/D. 
This  data  shows  considerable  scatter  at  the  lew  values  of  L/D* 
The  test  section  construction  procedure »  discussed  in  detail  in 

x3ndix  r,  was  such  tl  at  very  small  fillets  of  sclder  remained  at 
each  en'  ot  the  heated  length*     In  these  small  sections,  this  li- 
the accuracy  of  determining  heated  length.     It  can  be  seen  that 
approaches  a  niniHUM  value  as  L/r  is  increased,  and  that  CHP  is 
essentially  independent  of  L/     Cor  values  of  L/r  greater  than  hO* 
This  is  in  agreement  with  the  results  of  Bergles  (l,  11),     It  she. 
be  noteJ  that   ^he  value  of  CI3F  for  L/C  of  10  (the  generally  reported 
minimum  length  in  th     Russian  papers)  is  sere  1$%  greater  tl>an  the 
'  limura  value. 

Results  ot  the  investigation  of  Entrance  i  are  not  plotted* 

Holding  all  flow  variables  constant  and  t  ng  enly  ohe  form  oi  the 

entrance  tc   the  lieated  length  yielded  the   following  datai 
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Entrance  Form  Heat  Flux 

Sharp  iSdge  3.65  x  106  Btu/ft2hr 

1(5°  Chamfer  3.86  x  106  Btu/ft2hr 

Rounded  3.62  x  106  3tu/ft2hr 

Reentrant  3.55  x  106  Btu/ft2hr 

20- D  Calming  Length  3.82  x  106  Btu/ft2hr 

The  value  for  the  shr.rp  edge  is  a  mean  of  two  test  values,  the  others  are 
single  tests.     The  20-D  Calming  Length  was  the  standard  form,    used  for 
the  remainder  of  the  experimental  study,  the  value  listed  is  the  result 
of  several  tests  at  the  same  conditions.     Complete  data  for  these  tests 
are  lis'ed  under  runs  no.  33  through  no.  3?  in  the  listing  cf  data, 
Appendix  C. 

In  the  study  of  roaximum  heat  flux,  subccclings  on  the  order  of 
200  Btu/lbyrj  and  flow     rates  on  the  order  of  5  x  10°  lbm/ft2hr  were 
obtained  with  a  test  section  pressure  of  85  psia.     These  conditions  lead 
to  a  ("HF  of  almost  6  x  10^  Btu/ft^hr.      Exact  data  on  these  tests  are 
contained  in  the  listing  cf  data,  Appendix  C,  rune  through  numbers 
29,  30  and  31. 
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Discussion  of  Experimental  Results 

The  data  obtained   in  the  experimental  phas  udy  shows 

the   same  general  dependance  of  CHF  upon  the   other  variables  as  has 
been  reported  in  similar  studies.     This  data  is  perhaps  a  little  more 
well  ordered  than  that  normally  reported.     This,  of  course,   is 
heartening  to  the   investigators,   and  is  believed  to  be  the  result   of 
careful  experimental  procedures  and  constant  cognizance  of  the  flew 
stability  problem.     It  is  felt  that   instrument  and  reading  accuracies 
would   combine  tc  insure  an  experimental  accuracy  of  five  per-cent  or 
better. 

The  curves  of  CHF  verses  exit  subcocling  terminate  at  the  maximum 
subcooling  obtained.     This  m-xiraum  is  dependent  upon  the  test  section 
inlet  temperature  and  pressure.     The  curves  show  the  maximum  subcooling 
available  consistent  with  the  capacity  of  the  installed  heat  exchanger 
and  the  temperature  of  the  city  ccoling  water  (minimum  inlet  temperature 
"  h7°  F).     The  additional  l£  or  so  3tu/lbm  of  subcooling  which  could  have 
been  obtained  by  the  installation  of  a  complex  test  section  pre-cooler 
did  net  seem  tc  justify  the  additional  complication. 

The  curves  of  CHF  verses  exit  subcooling  show  a  minimum  value  of 
CHF,  signifying  a  change  of  flow  regime.     This  minimum  occurs  in  the 
subcooled  region  mainly  because  of  the  manner  of  determining  conditions 
at  the  pcint  of  CHF.     Fluid  conditions  were  determined  using  standard 
thermodynamic  variables,  which  are  valid  for  equalibrium  states,  but 
which  do  not  indicate  the  physical  conditions  with  inherently  non- 
equilibrium  Bubeooled,  forced  flow  boiling.     What  is  needed  is  a  better 
means  of  describing  the  physical  state,  one  which  would  include  some 
allowance  for  the  ncn-equilibH urn  conditions  as  they  actually  exist. 
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The  use  of  a  void-f rar  tion  variable  has  been  si 

(35) »  but  this  is  very  difficult  tc 

the  same  variables  as  CKF.     The  use  of  sub 

a  difference  in  enthalpy  or  te  re,  .be 

accompanied  by  a  realization  of  its  inherent  inability 

the  actual  physfeal  phenomena. 

At  the  very  hirh  values   o:   exit  subcoclirv 
aeeras  to  become  less  steep.     That  is  tc  say  that  the  .  of 

CRF  with  subcooling  decreases  with  increasing  subeoolir    .  rery 

slight  effect,  not  at  all  pronounced  at  low  system  pressures.     The 
high  pressure  data  of  Doroshchuck  sty*  Lar^tsman  (10)  and  Omatski  and 
Kiehigin  (5,  6)  shows  this  effect  rsore  clearly*     A  possible  explanation 
for  this  lies  in  analysis  of  the  effect  of  subcooling  upon  the  overall 
boiling  mechanism.     It  appears  that  once  high  values  of  subcooling 
are  reached',  and  a  large  recendensation  potential  is  established, 
further  increases  in  subcooling  have  less  effect. 

The  increase  of  CRF  with  decreasing  tube  diameter  is  evident  in 
figure  5.     It  has  been  suggested  that  void -fraction  increases  with 
decreasing  diameter.     This  is  due  to  the  manner  in  which  void-fraction 
Is  defined,  rather  than  rcme  phvsical  e;  ,     The  increase,  in  CHF 

Is  a  direct  result  of  the  requirement  for  continuity  of  mass  flow.     The 
size  of  a  vapor  bubble  is  dependent  mainly  upon  pressure.     In  small  tubes, 
the  bubble  occupies  a  proportionately  larger  part  of  the  crosn  section, 
requiring  that  the  liquid  phase  assume  a  higher  velocity  to  maintain 
continuity.     It  has  been  repeatedly  shown,  and  can  be  easily  explained, 
that  CHF  increases  with  increasing  velocity j  thus  the  increase  in  CKF 
with  decreasing  diameter. 

The  comparison  of  CHF  with  length-diameter  ratio  (figure  7), 
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shows  a  leveling  off  at  a  value  cf  L/  0. 

~s  enable  minimum  value-  tl 
where  one  is  attempting  to  insi 
variables  a  "e  separated. 

An  interesting  effect  of  changing  L/':  -  of 

CHF  verses  diameter  (figure  5)  The  curve        ant  tube  length 
re  .fleets  the  effect  of  simultaneously  chanr:     araeter  and  length* 
diameter  ratio*  At  the  larger  diameter*  the  data         nt  tub© 
length  reflects  a  substantial  decrease  in  l/D,  thus  the  greater 
Unfortunately,  the  limits  cf  the  apparatus  were  such  that  lai 
diameters  could  not  be  investigated  at  the  mass  flow  rate  selected  for 
this  study.  It  is  altogether  possible  that  the  carve  for  constant 
length  would  show  an  increase  in  CHF  for  larger  diameters;  i.e.  for  even 
greater  reductions  in  L/  . 

Comparison  of  CHF  data  for  two  tube  diameters  (figure  6),  shews 
that  diameter  is  one  of  the  r^cre  important  variables.  Within  the 
experimental  accuracy  of  this  study,  one  can  conclude  that  the  variation 
of  CHF  with  diameter  is  nearly  independent  of  subcoolir.  .  ,   least  in  the 
low  pressure  region.  This  requires  investigation  at  higher  system 
pressure,  where  the  effect  cf  pressure  becomes  more  pronounced,  and 
where  higher  levels  cf  subcocling  are  available. 

The  data  from  the  entrance  effect  study  shows  that  entrance  geometry  is 
unimportant  in  determining  CHF.  Here,  as  thrc     b  this  problem,  there 
is  a  high  degree  of  interaction  between  the  variables.  The  invest igation 
of  chat  1    i-ntrance  form  in  this  study  was  done  with  a  large  length - 
diameter  ratio  to  minimise  the  L/D  effect.  As  a  result,  variations  in 
flow  due  to  the  change  in  entrance  form  were  likely  damped  out  at  the 
point  of  tube  failure.  Decreasing  the  heated  1<:  ould  have  introduced 
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the  complication  of  the  L/~  effect*  safe  tc  conclude  that  changing 

et  for  values  of  L/D«     The  somewhat  lower 

form  n  bility  due  tc  compressibility 

..■be  non-flow  spare  (  tionf  but  is  mere  lik< 

to  expt.  /oer, 

.  ff»heat  flux  data  is  ci  interest  only  for  future  investigations 
ueing  the  apparatus  employed  in  this  study*     This  has  no 

.he  overall  CHF  problem  because  of  the  limitations 
imposed  upon  pressure  rate  by  the  apparatus.     Tests  were 

to  determine  the  maximum  CHF  that  could  be  cbtaj  Lth  vjorking 

size  tube  disaster  (0»2l2  in.)  and  this  apparatus. 
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The  following  cowmen  rioua  cor  ieting 

schemes  are  based  upon  close  examination  air  ,h  experii 

data.  The  effect  of  the  controlling  parameters 
described  in  detail  in  the  section  in  Factors  Influei 
Heat  Flux  (Chapter  III). 
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3ellt     (ref;   21) 

The  use  oi  this  cci"  restricted  t  sia. 

There   is   no  restriction  on  the  tube  diameter,  but  fr  r  results, 

it  has  been  shown  C  eterg 

therefore,  this  correlation   should  be  considered  useful  for  lsa 
diameters  (0.25  ^in.    or  larger)  only.     Although  the  accuracy  is  \$%  to 
in  the  case  illustrated   (figure  9),  the  e  fc,ion  is  relatively 

simple  and  easy  to  apr-ly,  and  is  generally  conservative  in  it  diction* 

Jens  f_  Lottes;     (ref:  17,) 

This  correlation  is  of  a  power  function  form,  and  would  indi- 
zero  CHF  at  zero  Bubcooling  unless  a  realistic  lower  limit  was 
established,     ^ince  the  correlation  was  des  to  have  no  negative 

deviation,   it  i  sf  the  most  conservative  consider*  main 

virtue  is  its  ea^e  of  application. 
Ryai  oy  and  Berzina ;     (ref;  36) 

The  2  vm  correlation  is  adequate  for  the  data  at  pressures  above 
lh7  psia.     It  is  not  consistent  with  ex  nt  for  the  lower  ■  re 

data  (30  or  90  psia).     There    ■  Imm  subcocling,  approximately 

18  Btu/lbrc  below  which  the  correlation  fails  due  to  the  r  si 

form. 

The  correlation  for  B  *  hmrr  is  adequate  for  larger  tubes, 
(approximately  .3 in. ),  but  for  smaller  tubes  it  suffers  from  lack  of 
consideration  i  iaraeter  effect*     Ryabcv  and  Berzina  considered  the 

diameter  effect  in  the  derivation  of  this   correlation  by  changing  the 
leading  coefficient,  but  they  intimated  that  the  diameter  effect  is 
negligible  for  larger  diameters.     This  correlation  could  be  improved 
by  substituting  a  function  of  diameter  for  the  leading  constant. 


3U 


Zenkevicht      (reft  ih) 

viscosity  in     his  coi  ture 

produce?  a  curve  scmowh&t  li  nded 

if   this  correlation  is  tc  be  used, 

ce  cf  a  tube  -Hamster  factor  limits  the  us^  af  this 

correlation.     7ne  lower  limit  of  0,16  in.   is  toe  small  tc  exr-iude  all 
diameter  effect.     It  is  recc  id  that  the  use  of  this  cori 

be  limited  to  tubes  of  0,2: '<  in.   in  diameter  cr  larger,  Where 
dependence  of  CHF  u  Laraeter  Is  quite  sr  all. 

Zenkevich  &  Subbotin;     (ref?  Ik) 

This  correlationf  lil  neglects  the  effect  of  tube 

~ter.     The  lower  limit  of  0.16  in.   is  too  small  to  neglect  this 
effect  entirely.     For  the  larger  diameter  (.316  in.)  data;   this 
cct  on  is  excellent.     The  lower  limitation  on  subeooiing  is  sound, 

■due  to  the  vv  -.Jcal  forsra  of  the  equation.     The  upper  limit  has 

nd  seess  evidence  of  caution.     In  spite 

cf  the  diameter  consideration,  the  mathe  1  expression  is 

convenient,   a  S  results  f c  -  ;-r  tubes  are  quite  good. 

Gambill:     (re  is  3,0) 

The  curve  for  wall  t  ure  is  suspect,  and  available  data  of 

wall  superheat  at  burnout  lees  not  bear  out  the  graph  presented.     There 
is  too  much  latitude  in  choice  of  constants  and  exponents  to  be  used 
in  the  correlation.     Presumably,   if  the  ■  cnablv  well  ordered, 

the  fit  of  the  correlation  can  be  improved  by  adjusting  the  values  of 
K,  K',  n,  and  n.     There  are  no  easily  applied  guidelines   for  the 
adjustment  ©I  these  constants  included  in  the  presentation,     ft 
correlations  does  include  8  ter  ter»,  but  it  is  only  in  the  forced 

convention  part,  and  then  to  onlv  the   (m  -  l)  power,  where  m  is   one  of 
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Gambill's  "variable"  c ens  ants ,  generally  on  the  order  of  0.8,  This  is 
one  ol  the  few  correlations  which  includes  small  tube  diameters, 
clearly  the  correlated  effect  of  diameter  is  loss  than  that  observed  in 
the  data  of  Kergles  (l),  Poroshchuck  and  Lantsman  (9),  and  figure  lw 
The  correlation  is  no  better  than  that  of  Mirshak  (13)  at  lew  pressures 
and  is  considerably  clumsier  to  handle.  At  high  pressures  (figure  9) 
this  correlation  is  no  better  or  worse  than  many  ethers. 
Be  math:  (ref:  16) 

For  a  riven  pressure,  tube  diameter  and  flew  rate,  this  correlation 
is  linear  with  subcooling.  This  is  a  faitly  good  first  approximation 
and  typical  cf  many  correlations.  This  correlation  is  rather  difficult 
tc  apply,  requiring  auxiliary  calculations  tc  determine  heat  transfer 
coefficients  and  wall  temperatures.  Of  the  seven  high  pressure 
correlations  plotted  en  figure  9$   this  is  one  of  the  best,  predicting 
CEF  within  $%   for  the  selected  conditions, 
3ettis»  (ref ?  17) 

■  —■■*■■■———      n  n  —   ii  wii  ■     ■  ■  ■  Hi 

The  use  of  this  correlation  is  restricted  to  pressures  near 
2000  psia.  There  is  no  restriction  on  a  minimum  diameter;  therefore, 
the  corr  laticn  will  yield  low  predicted  values  of  CRT  for  diameters 
less  than  .25  in. 
Griffith t  (ref:  13,  17) 

This  correlation  was  ogiginally  determined  for  several  fluids 
including  water  and  it  is  applicable  throughout  the  subcooled  boiling 
region  and  into  net  steam  generation.  Generally,  over  a  wide  range 
of  variables,  this  correlation  is  the  most  consistently  accurate  one 
available.  Its  disadvantage  lies  in  its  extremely  cumber.0 c.^e  and 
unwieldy  form.  Computation  can  be  simplified  screwhat  by  using  the 
curves  presented  in  ,iIAFD-.lS8  (17)  for  several  cf  the  recurring  variable  groups. 

36 


T-iirshak,  et  al:      (refi  13) 

The  form,  cf  the  correlation  is  linear  ■ 
exhibits  the  approximate   to  ental  da 

the  range  of  :he  corre] 

dependence  of  CHt    en  diameter  i:  cted.     The  follow:  neral 

cements  are  appropriate: 

1)  The  correlation  rives  a  conservative  indication 
of  CHF,  but  errs  by  7$%  - 

2)  There  is  no  consideration  of  the  variatic 
CPF'  with  diameter 

3^     The  correlation  is  safe,  easy  to  apply 

■nt:     (^ft   6) 

This  correlation  was  designed  for  correl  ba  for 

tubes  of  a  specific  length.     Hew  te  Subcoolin?  Gradient  for 

ether  tube  lengths  is  unresolved.      Subecoling  gradient  is  not  an 
independent  variable  equation,   b   Ln      bied  to   Heat  Flux  by  the 

First  taw  of  Thermodynamic?.,   Diameter,   and  flow  rate.     No  diameter 
effect  is  noted  between  0.30  -  0.88  in.     This  is  a  reasonable 
approximation.     Variation  in  pressure  between  60  -1000  psia  is 
considerable  and  should  be  taken  into  account.      Since  the  equation 
was  derived  from  e  LI  sample  of  data  and  for  tubes  of  only  one 

length,  its  use  is   net  recommended. 
Doroshchack  <^  Lantsman:      (ref;  9) 

This  is  one  of  the  few  studies  which  has  been  devoted  exclusively 
to  the  diameter  effect  problem.     The  general  trend  is  correct,  but 
comparing  this  correlating  equation  with  experimental  data  -  even  that 
from  the  yaper  in  which  this  correlation  Is  ^resented,  shows  an  error  of 
up  to  30/'.     Part  of  the  difficulty  is  due  to  the  interaction  of  the  variables. 
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s  soh<  bably  that 

in  (l,   5,   ' ,)  -  s  (1,  11,)  owed 

the  cc^  ceed  a  cror.  w  of 

A   [(d/Dq)1-2  -  lj     fit  the  OmatsJd  and  Kichigin  d:  ,   $  and  6) 

quite  well,  but  this  failed  in  comparison  with  other  studies. 

It  is  net  rtsccrrriended  that  this  cor  'ct 

j  however,  it  will  serve  tc   apply  a  diameter  eration  to  ct 

correlations,  -within  rather  bread  limits. 
KTT  Srall  ri arret er 

This  equation  was  developed   not  sc  s  a  prediction  of  CHFf  but 

tc  provide  a  fairly  concise  mathematical   expression  for  ti  pressure 

lata.  *Tc  pressure  term  is  ir?  laded  as  the  variation  of  CHF  with 
pressure  in  neither  well  ordered  nor  of  large  magnitude.  The  change 
in  TFT  with  pressure  over  this  range  '       i  than  the  error  of 

tt'r  icn  due  to  1  sumption   of  '  e   linear  function  of 

sub.-    sli      .  '  s,  in  fact,   not  li*  ■■  ith  subcooling,  this   form 

has  been  assumed  merely  for  convenience.     It  is  he 

use  of  ^his  equation  be  limited   tc  dr  on  of  general  levels, 

and   that  the  actual  data  be  used  in  any  attempts  at  accurate  prediction  of 
CHF. 


•IX  I 

PR0C 

Figures  3,  ii,  10,  11  and  12  provide  a  graphical  presentation  of 
CKF  for  pressures  from  30  psia  to  2&U0  psia  ever  a  range  of  mass  flow 
rates  and  exit  s^becolings  for  tubes  of  two  ram  diameter.     These  curves 
in  conjunction  with  the  data  of  CKF   dependence  upon  tube  diameter 
provide  a  means  of  graphically  determining  CHF  for  tubes  up  tc  one- 
fourth  inch  diameter,  flow  rates  from  2  x  106  lb^ft2  hr.  to  22  ItWft2  hr, 
at  any  pressure  and  exit  subocx  ling.     A  procedure  of  linear  and  graphical 
interpolation  is  recommended   lor  predicting  CHF  for  any  reasonable 
small  diameter  application.     In  eases  where  linear  interpolation  is 
recommended   ,  auxiliary  curves,   have  shown  this  technique  tc  be 
satisfactory.     The  accuracy  of  this  graphical  prediction  technique 
will  exceed  that  of  the  few  correlating  equations  applicable  to  email 
tubes,   and  in  the  majority  cf  the  eases  tested,  was  within  13%  of  the 
experimental  data  considered. 

The  following  procedure  is  recommended  for  predicting  CHF: 

1.  Select  desired  system  pressure,  mass  flow 
rate,  tube  diameter  and  length. 

2.  Locate  the  plots  for  pressures  bracketing 
the  desired  pressure. 

3.  At  the  selected  exit  subcooling,  read 
CHF  for  the  three  or  four  mass  flow  rates 
presented,  for  each  pressure. 

it.     Prepare  interpolation  curves  with  coordinates 
of  CHF  and  Mass  Flow  Rate,  and  determj  ne  CHF 
corresponding  to  the  desired  flow  rate  at  each  pressure. 
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$.     Using  the  two  values  thus  obtained,  Unear 

interpolation  between  the  pressures  u 

consideration  will  give  a  value  of  CEF  for 

the  desired  pressure,    flow  rate  and 

subrooling,   for  a  diameter  cf  two  mm. 

Refer  to  figure  5,  ising  the  curve  for 

constant  L/D.   determine  the  value  of  CPF 

for  the    ;iameter  in  question.     The  ratio 

of  this  last  value,  to  that  at  two  ran? 

(a  value  of  CHF  of  li.l£  x  106  Btu/ft2  hr. 

at  2  mm  diameter  will  provide  satisfactory 

results)  multiplied' by  the  previously 

determined  value   cf  CHF,  will  yield  a 

predicted  value  of  CHF  at  the  desired 
conditions  cf  pressure,  flow  rate,  subcocling 
and  tube  diameter. 

7.     In  cases  where  the  desired  L/D  is  less  than 
20,  the  same  sort  of  ratio  correction  should 
be  performed,  using  figure  7.     For  values  of 
I/D  of  20  or  greater,  the  error  in  CHF  will 
be  tec  small  to  warrant  consideration. 
For  a  rough,  quick  estimate  of  Cffi',  visual  interpolation  for  mass 
flow  rate  on  the  plot  for  pressure  nearest  the  desired  pressure  can  be 
easily  accomplished.     The  diameter  correction  should  be  included  in  this 
abbreviated  procedure.     This  same  sort  of  approach  is  useful  in 
estimating  the  effect  of  gross   change  in  the  important  variables. 
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Experimental  on  was  conducted  en  exis  which 

been  constructed  for  a  previot r  been  used 

fcr  several  other  boilirij  .  r  miner 

modification  was  neessary  for  this  inve  -ion,      The  apparatus 

is  shown  schematically  in  figure:  1.     The  fluid  loop  consists  of  a 
closed  main  circulation  loop,  with  a  by-pai  rallel  test  section 

line.     All  components  are  of  corrosion  resistant  materials, 
system  contains  a  main  circulating  pump,  accumulator,  the  test  section 
line  with  flowmeter,  preheater,  test  section  with  instrumentation, 
main  loop  by-pass  line  and  heat  exchanger  with  city  water  cocli     . 
Auxiliary  to  this  are  a  fill  pump,  supply  tank,  degassi 
continuously  lemineraliter*     Power  is  supplied  to  the  test 

section  by  a  pair  of  a.c.  motor  -  d.c  .         orator  sets* 

The  main  supply  pump  is  a  two-stage  turbine  type  pump,  driven  by 
a  3  hp  induction  motor*     The  purp  is  ca;  ivering  3*6  gallons 

per  minute  at  260  psi.     A  blad  pa  arrumulatcr  locate-  he 

exit  side  of  the  pump  serves  to  '-amp  out  pressure  fluctuations*     Plug 
valves  at  each  end  of  the  by-pass  lit  -cntrcl  for  the  test 

section  pr         re*     Test  section  flow  rate  is  determined  by  a  Fisher 
Pointer  flowmeter  '  rd  controlled  with  a  Hoke  needle  valve.     This  needle 
valve  also  provides  the  15>0  psi  up-stream  ire  drop  necr  be 

insure  flow  stability. 

The  mctcr-generr.tor  sets  consist  of  U<0  v.,  3  phase  synchronous 
motors  and  36  kilowatt.  d,c*  generators,  connected  in  series*     lis.eh 
generator  is  nomijially  rated  at  12  volt  and  3000  amps.     A  total  power 
of  approximately  60  kilowatt  is  available  at  the  test  section* 
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rumentation  is  available   i         onitoring  pressure  levels, 
temperatures j  test  section  flow  rate  and  test  section  current  and   voltage. 
Test  section  exit  pressures  are  read  on  Bourdon  type  pressure  gauges 
with  specified  accuracy  of  0.5>  psi  or  better*     The  gauges  are  calibrated 
periodically  with  a  dead  weight  tester.     Test  section  flew  rate  is 
determined  with  t^e  Flowrater  "etersj  the  range  of  flow  rate  is  varied 
by  changing  tubes  an     floats,      ''iters  are  calibrated  in  place  and 
calibration  curves  are  available.     Test  section  power  is  obtained  from 
the  test  section  voltage  and  current.     Voltage  across  the  test  section 
is  read  directly  with  a    feston  variable  scale,  calibrate:]  d,c,  voltmeter. 
Current  is  read  across  a  national  Bureau  of  Standards  shunt,  with  a 
specified  calibration  of  60,17  amps/mv. 

All  temperatures  are  determined  by  thermocouples  manufactured 
frc~  30-gauge,  duplex,   Ccpper-Ccnstantan  thermocouple  wire.     Test 
section  inlet  and  exit  temperatures  are  measured  on  thermocouples 
inserted  into  the  flow  v.ith  Conax  fittings.     Output  voltages  of  the 
thermocouples  and  current  shunt  are  displayed  upon  a  Minneapolis* 
Honneyv all      1'cwn,  0-26  mv,,  single  channel  recorder, 

A  much  more  detailed  description  of  the  apparatus  is  contained  in 
ref,   1}  however,  certain  modifications  have  been  made.     The  mere  important 
of  these  include: 

a.     The  braided  copper  flexible  power  connection 
anrJ  the  rigid  up-stream  jumper  iron;  the 
test  panel  to  test  section  bus  connection 
were  replaced  during  this  study  by  liaekwerth 
Eeese  Flexible  Dry  Jumpers,  This  requires 
that  the  test  section  be  supported.     Support 
is  pcrvided  by  a  laboratory  support  and  clarnp, 

61 


with  electric al 

b.     The  braided  copper  y  cables  have 

beenrepLiK  iese,  water 

coolod  cables*     These  cal  re  cooled  by 

city  water  tapped  off  of  the  he 
coolant  supply  lin  1  to  a  drain, 

e.     The  preheater  haw  bee-  y  the 

addition  of  three  $  kilowatt  iusnersion  heaters, 
with  separate  controls.     These  heaters  are  in 
sertaes  with  the  original  6  kilowatt  variable 
control  heater  and  de  any  desired 

preheating  from  0  to  21  kilowatt,  with 
steady  state  operation. 

d.  The  steam  heater  in  the  suppl- 
been  removed  as  it  was  found  that 
sufficient  heating  and  debase '  uld  be 
obtained  froir  the                           eaters  alone, 

e.  The  thermocouple  selection  switch  now 
consists  of  two  Leeds  and  Horthrup,  12 
position  selector  switches,  in  series* 
Thermocouple  wires  are5  connected  through 
screw  type  t  1  boards,  thus  reduci 

t-up  tiro,     Th:'         "  fie  at  i  en  was  raide 
in  order  fee  iaerease  the  number  of 
thermocouples  which  could  be  patched 
to  tfas  ->.     Bcr ■■"         ee  tests  sho. 

ao  apprv  as  a 

result  of  the  added  .     In  this  study, 
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only  the  primary  srwj 

f.     When  usin  hort  test  sr 

or  those  of  type  el,   the  test 

section  re 

nfikc  ■  enerator 

highly  erratic.     To  alleviate  this 

blew,  a  series  resistance  w 
fabricated.     This  resistance  sted 

of  a  10  inch   length  of  0.2L2  inch  t 

Lnless  steel  tubing ,     '   h  bus 
connections,  place-  electrically  in 
series  with  the  tc  coded 

v,  r  „  -H-.  -  eater*    This  a 
t      sed  an  e 

apprcscir  "  .s  across  the 

supply,  -which  inprcved  the  r..»r .    set 

Tirol  considerably,     Towe- 
section  ed  by  readiis 

drop  across  the  test  section  only,  and  by 
reading  current  th  ries 

'  en. 

f.     In  the  tests  where  very  high  power 

applications  -wore  required  tc  reach  CHF, 

the  tost  section  exit  bus  connection  was 

found  to  heat  op  excessively,  causing 

■ptokii  in  t        a  snrcal!      Li • "  ,   not 

tc  n-  • 

test  sections.  Tn  *ien  cf  this,  a  bus 

connection  coolant  scher^e  was  devised.  A 
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-    ich  copper  nipple  and  two  elbows  were 

silver  brazed  tc.  the  brass  bus  connection, 
water,   fro?  the  power  lead   pooling 

system,  was  diverted  through  this  nipp     , 

drawing  heat  from  the  bus  connection.      Heat 

balance  checks  when  operating,  comparing 

power  input  to  fluid  temperature  rise  through 

the  test  section  showed  that  this  cooling 

device  had  no  effect  on  the  test  data.     No 

further  overheating  problems  were  encountered 

once  this  cooler  was  installed,  and  it  was 

used  whenever  there  was  a  requirement  for 

high  power  input. 
Test  Section 

Because  each  test  destroyed  the  test  section,  construction  was  kept 
as  simple  as  possible.     A  test  section  was  constructed  by  facing  a  3 A  inch 
brass  pipe  cap  and  drilling  a  hole  in  the  center  of  it.     A  section  of 
the  test  tubing  was  silver  soldered  into  this  hole.     At  the  other  end 
of  the  desired  heated  length,  a  5>/l6  inch  brass  bushing,  about  £/8  inch 
long,  drilled  to  the  O.D  of  the  tubing  was  silver  soldered  to  the 
tubing.     Heated  length  was  measured  from  fillet   bo  fillet,     A  calming 
length  of  approximately  20  diameters  extended  beyond  the  up-stream 
bus  bin.  . 

Sections  were  installed  in  the  apparatus  by  threading  the  pipe  cap, 
the  down-stream  end,  to  the  down-stream  plenum  chamber  which  contained  the 
exit  pressure  tap  and  exit  thermocouple.     The  up-stream  end  was  installed 
by  inserting  the  calming  length  into  a  Conax  fitting  sealed  with  a 
neoprene  washer.     Brass  bus  connections  were  installed  on  the  bushing  and 
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In  the  st  bine  the  out 

diameter  o  ,  obtain 

adequate  ::  me  ros  .  -  noted  in  the 

,ing  or  'Hate  ,  runs     ^9,  ':■ 

Be  ryll  iu^-C  oppe  r  Tu 

At  the  beri  -f~u 

was  order-:    ,  .  ra)«     It  "elt 

that  this  tu1'  stainless 

stee  .  »st  sections 

were  wade  up  ere  run,  '.ffir-uity 

ecuaterod  using  the  Be-O  .  it 

problem  a*a  . 

About  the  only  fabric  tubing  hi 

that  It  could  be  e~~  .     Unfortunately,  the 

alloy  was  not  B  D  degress  F. 

test  sections  \r  soldering,  tJ  near 

-iclder  tions  underwent  a  phase  t  nied 

by  a  el  in  electrical  resir  I      .    Since  this  change 

in  re;:-  red  only  in  the  region  where  heat  vs&B  applied  for 

solder?-   ,  -used  as  uneven 

distri  .  ~t 

section  voltage  drop  and  current  did  net  :' 

Test  sections  were  solder,  but  the  h: 

exit  t  •    /ares  enc  "d  the  soft  solder 

to  puddle  j  uat  '.  .  leaks  las 

iel+  less  -    .  -ature  -  oft 

Mer  was  tried,  with  r. .  r   results. 

65 


■ 


The  lower  resistance  Be-Cu  tab:  ised  m.g.   set  ccntrcl  probL 

The  series  resistance  previously  mentioned  alleviated  this  problem,  but 
tended  to  oornpl:'  cate  the  procedure. 

data  was  obtained  using  this  Be-Cu  tubi-    ,     Subsequent  comparison 
with  sirilar  data  obtained  with  type  30ii  stainless  steel  tubing  showed 
that  the  Be-Cu  data  was  inconsistent*  for  either  the  soit  soldered  or 
silver  soldered  sections. 

A  corrosion  check  ves  made  with  the  Be-Cu  tubing.     After  one  hour 
at  a  heat  flux  of  a  I  cut  one-third  of  CKr ,  the  inside  of  the  tube  showed 
distinct  blackeninr,   using  derassed  and  demineralized  water. 

In  view  of  the  various  difficulties,  it  was  concluded  that  the 
Be-Tu  tubing  was  unsatisfactory  for  experiments  o;  this  type,  and  thick 
wall,  type  30b  sr,ainlesn  steel  tubing,  with  inside  diameter  of  0.079 
inches  (2  bib)  was  obtained. 

The  stainless  steel  tubing  could  be  silver  soldered  satisfactorily 
and  the  data  obtained  iror  it  shewed  considerably  rcre  consistency. 
The  data  presented  as  the  result  of  this  investlgs  ion  was  obtained  using 
type  30L  stainl  ■     stoel,  except  where  noted  to  be  Be-ru  (runs  02  through 
#23 )l  and  except  for  the  D  ■  0.180  inch  tubing  which  was  type  "A"  Nickel. 
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Often  the  process  of  obtaining  the  desired  exit  subcooling  was  a 
very  delicate  one.  This  was  particulary  true  whe 

versus  subcooling  was  approximately  parallel  to  the  curve  of  heat  flux 
versus  subcooling  as  determined        irst  Law  of  Thermodynamics. 

In  the  large  majority  of  the  test  runs,  it  was  possible  to  obtain 
the  desired  valua.  of  exit  subcooling  within  $%   at  the  point  of  CHF. 
In  these  few  cases  where  exit  subcooling  varied  more  than  $%   from  the 
desired  value,  an  auxiliary  plot  of  CHF  versus  subcooling  was  made,  to 
deterrd-me  the  appropriate  value,  or  to  determine  that  the  variation  of 
CHF  with  subcooling  was  not  significant  in  the  range  of  subcooling  in 
question. 

For  the  purpose  of  this  study,  CHF  was  defined  as  the  heat  flux  at 
tube  failure.  While  the  heat  flux  thus  determined  may  exceed  that  at 
the  departure  from  nucleate  boiling  (DNB),  heat  flux  sc  defined  should 
coincide  closely  with  the  hest  flux  at  the  point  of  transition  from 
nucleate  to  film  boiling.  Mo  burnout  protection  device  was  used  in 
this  study,  as  it  was  felt  that  less  time  and  jreterials  were  involved 
in  the  fabrication  of  a  new   test  section  for  each  run  than  in  the 
construction  of  such  a  device.  Furthermore,  the^e  can  be  no  doubt 
that  a  genuine  maximum  heat  flux  has  been  established  when  the  tube 
ruptures. 

All  experiments  were  conducted  with  the  test  section  in  the 
horizontal  position.  System  water  was  degassed  only  periodically. 
Care  was  taken  to  insure  that  no  air  was  trapped  in  the  oystem  when 
replacing  sections  and  extensive  air  bleeding  was  done  before  commencing 
a  new  run.  The  inside  surface  of  the  test  section  was  given  a  thorough 
cleaning  with  acetone  prior  to  installation.  No  additives  were 
introduced  into  the  water. 
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In  almost  all  cases,  the  tube  ruptured  within  three  diameters 
of  the  exit  connection.     Thus  renditions  in  the  exit  plenum  chamber  were 
considered  to  be  representative  .  the  point  of  failure.     A 

recent  study  by  Lopina  (3)  has  shown  that  the  pressure  within  the  plenum 
is  virtually  identical  with  that  at  the   down-stream  end  of  the  heated 
section  throughout  the  subcooled  region.     The  test  section  and  plenum 
were  wrapped  with  fiberglass  insulation  to  greater  than  the  critical  radius 
to  minimize  heat  lot    . 

Because  the  recorder  used  presents  only  one  channel,   it  was 
necessary  to  estimate  certain  quantities  at  the  point  of  CHF;   however, 
near  failure,  very  small  increments  of  power  were  a  ,  and  all 

quantities  were  checked  frequently.     Data  not  obtained  right  at  failure 
could  be  readily  inferred   from  that  recorded  just  prior  to  reaching 
CHF.     All  changes  in  conditions  were  made  very  gradually  during  a   run, 
and  once  the  procedure  was  well  established,   it  was  possible  to  obtain 
a  set  of  data,  e.g.  one  data  point,  in  approximately  forty-five  minutes. 
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APP^TOIX  F  -  SAXPLE  DATA  SHEET 

D  O.Q1S>      in.        P        5Q  psia     Run  #         4S 

L  3.15        ln.     L/D       4Q  Date  ZflS 
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CKF  =      1C?.55-k12^S»(  60.17)  (5.41  )        =     4.BZ   BtuA.,2 

Obtain   Temperature   from   Thermocouple   Conversion  Chart 
h  sat  =   h  at  p  h  bulk  =   h  at  Tout 
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APPENDIX  G 
DEVELOPMENT  OF  THE  MIT  SMALL  DIAMETER  CORRELATION 

This  equation  was  developed  by  the  authors  for  the  low  pressure 
MIT  data  and  fits  the  majority  of  it  to  *  2$%, 

(q/A)cr  «  1.738  x  106  (U.U8  +  .Ol^hsat)  W%)***  (D/D0)-«^/i/d\-.12 

The  correlation  was  developed  in  the  following  manner:     The  30 
and  90  psia  plots  (figures  3  &  U)  for  2  mm  data  (.079  in.)  were  linear- 
ized to  determine  the  slope  and  the  intercept  of   the  part  of  the  equation 
dependent  on  subcoolong.     The  flow  rate  dependence  was  determined  by- 
plotting  on  log-log  paper,  cross-plots  Of  the  2  ran  data  at  constant 
subcoolings.     Also  plotted  were  cross-plots  of  Bergles!(ll)  and 
Vfessel's  (26)  data  for  constant  subcoolings.     A  mean  value  of  the  slopes 
was  used  to  determine  the  exponent  of  the  flow  rate  term.     The  same 
technique  was  used  for  diameter  and  L/D  dependence,  using  the  data  of 
figures  $  and  7,  and  the  data  of  Bergles  (11).     The  leading  constant 
was  determined  by  the  zero  subcocling  intercepts  and  the  values  of 
the  constants  G0  and  D0#     No  pressure  term  is  included  as  the  variation 
of  CHF  with  pressure  over  -this  range  is  neither  well  ordered  nor  of 
large  magnitude. 

The  correlation  was  plotted  over  the  MIT  data  and  small  adjustments 
were  made  in  the  constants  and  exponents  to  get  a  better  all  around  fit 
and  balance  out  the  positive  and  negative  deviations. 

A    f  urther  refinement  of  the  correlation  would  be  to  make  the 
heat  flux  a  weak  2nd  or  3rd  order  dependence  on  subcooling,  eliminating 
the  linear  approximation.     A  minimum  could  be  introduced  at  about 
Absat  •  20  to  30  Btu/lbm. 
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APPENDIX  H  -  NOMENCLATURE 

—■■«»—■■     i  ii ■    i wii.i    imui«niw>  ■ — ■■h»i.i— «»— — — 

t  -  Acceleration  in  g*s 

p 

A  -  Heat  Transfer  Area,  ft 

c  -  Specific  Beat,  Btu/lbm-deg  F 

D  -  Tube  Diameter,  inches 

g  -  Acceleration  of  gravity,  ft /sec2 

go  -  32.2  lbf-ft/lbm-sec2 

h  -  Heat  Transfer  Coefficient,  Btu/hr-ft2-deg  F 

h  -  Enthalpy,  Btu/lbm 

J  -  Mechanical  Equivalent  of  Heat,  778  ft-lbf/Btu 

k  -  Themal  Conductivity,  Btu/hr-ft-deg  F 

L  -  Tube  Length,  inches 

p  -  Pressure,  psia 

Pr  -  Prandtl  Number 

q/A  -  Ifeat  Flux,  Btu/hr-ft2 

t  -  Temperature,  deg  F 

T  -  Teaperature,  degrees  Eankine 

ATsat  "  (Tsat  "  Tb)»  &*&  F*  Subcooling 

v  -  Specific  Volone,  ft^/lbm 

V  -  Velocity,  ft /sec 

**  -  Viscosity,  ft2/hr 

M'  -  Viscosity,  Ibm/hr-ft 

(*  -  Density,  lbm/ft3 

*""  -  Surface  Tension,  lbf/ft 

Ah  •  (hgat  -  life),  Btu/lbm,  Subcooling 
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Subscripts 

bo  -  At  the  point  of  burnout  or  critical  beat  flux 

CT  -  At  the  point  of  burnout  or  critical  heat  flux 

b  -  Bulk  or  raided  isean  condition 

crit  -  Critical* 

fg  -  Saturated  Condition  •  Vapor  Condition 

g  -  Vapor  Properties 

1  -  liquid  Properties 

pb  -  Pool  Boiling 

•  -  Saturated  Conditions 

sat  •  fiaturafced  Conditions 

▼  -  Vapor  I'roperties 
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